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Abstract. - The Landau Zener method allows to measure very small tunnel splittings A in 
molecular clusters Fes ■ The observed oscillations of A as a function of the magnetic field applied 
along the hard anisotropy axis are explained in terms of topological quantum interference of 
two tunnel paths of opposite windings. Studies of the temperature dependence of the Landau 
Zener transition rate P gives access to the topological quantum interference between excited 
spin levels. The influence of nuclear spins is demonstrated by comparing P of the standard 
Feg sample with two isotopically substituted samples. The need of a generalised Landau Zener 
transition rate theory is shown. 



During the last few decades, a large effort has been spent to understand the detailed dynamics 
of quantum systems that are exposed to time-dependent external helds and dissipative effects 
[Q. It has been shown that molecular magnets offer an unique opportunity to explore the 
quantum dynamics of a large but finite spin. These molecules are the final point in the 
series of smaller and smaller units from bulk magnets to single magnetic moments. They 
are regularly assembled in large crystals where often all molecules have the same orientation. 
Hence, macroscopic measurements can give direct access to single molecule properties. 

The most prominent examples are a dodecanuclear mixed-valence manganese-oxo cluster 
with acetate ligands, Mni2 0, and an octanuclear iron(III) oxo- hydroxo cluster of formula 
[Fe8C>2(OH)i2(tacn)6] 8+ , Fes [pf, where tacn is a macrocyclic ligand. Both systems have a spin 
ground state of S = 10, and an Ising-type magneto-crystalline anisotropy, which stabilises the 
spin states with the quantum numbers M = ±10 and generates an energy barrier for the 
reversal of the magnetisation of about 67 K for Mni2 and 25 K for Fes PL [§■ 

Fes is particular interesting for studies of quantum tunnelling because it shows a pure 
quantum regime, i.e. below 360 mK the relaxation is purely due to quantum tunnelling, and 
not to thermal activation El. We showed recently that the Landau Zener method can be 
used to measure the very small tunnel splittings A in Fes M- The observed oscillations of 
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A as a function of the magnetic field applied along the hard anisotropy axis are explained in 
terms of topological quantum interference of two tunnel paths of opposite windings which was 
predicted by Garg ||] . This observation was the first direct evidence of the topological part of 
the quantum spin phase (Berry or Haldane phase H, [§]) in a magnetic system. 

Recently, we demonstrate the influence of nuclear spins, proposed by Prokof'ev and 
Stamp ||, by comparing relaxation and hole digging measurements ]Io| ] of two isotopically 
substituted samples: (i) the hyperfine coupling was increased by the substitution of 56 Fe 
with 57 Fe, and (ii) decreased by the substitution of 1 H with 2 H. These measurements were 
supported quantitatively by numerical simulations taking into account the altered hyperfine 
coupling JXy, 11 . 

In this letter, we present studies of the temperature dependence of the Landau Zener 
transition rate P yielding a deeper insight into the spin dynamics of the Fes cluster. By 
comparing the three isotopic samples we confirm the influence of nuclear spins on the tunneling 
mechanism and in particular on the lifetime of the first excited states. Our measurements 
show the need of a generalised Landau Zener transition rate theory taking into account 
environmental effects such as hyperfine and spin-phonon coupling Jl^| . 

All measurements of this article were performed using a new technique of micro-SQUIDs 
where the sample is directly coupled with an array of micro-SQUIDs | fl3"l| . The high sensitivity 
of this magnetometer allows us to study single Fes crystals jl4| of the order of 10 to 500 /im. 

The crystals of the standard Fe8 cluster, st Fes or Fes, [Fes(tacn)602(OH)i2]Brs.9H20 where 
tacn = 1,4,7- triazacyclononane, were prepared as reported by Wieghardt et al. H|. For the 
synthesis of the 57 Fe-enriched sample, 57 Fes, a 13 mg foil of 95% enriched 57 Fe was dissolved 
in a few drops of HCI/HNO3 (3 : 1) and the resulting solution was used as the iron source in 
the standard procedure. The 2 H-enriched Fes sample, D Fes, was crystallised from pyridine-ds 
and D 2 (99%) under an inert atmosphere at 5°C by using a non-deuterated Fe(tacn)Cl3 
precursor. The amount of isotope exchange was not quantitatively evaluated, but it can be 
reasonably assumed that the H atoms of H2O and of the bridging OH groups, as well as a part 
of those of the NH groups of the tacn ligands are replaced by deuterium while the aliphatic 
hydrogens are essentially not affected. The crystalline materials were carefully checked by 
elemental analysis and single-crystal X-ray diffraction. 

The simplest model describing the spin system of Fes molecular clusters has the following 
Hamiltonian 

H = DS 2 2 +E(S 2 x -S 2 y ) +H 2 -gn B n Q SH (1) 

S x , S y , and S z are the three components of the spin operator, D and E are the anisotropy 
constants, H 2 takes into account weak higher order terms Jl5|, |l(|, and the last term of 
the Hamiltonian describes the Zeeman energy associated with an applied field H. This 
Hamiltonian defines a hard, medium, and easy axes of magnetisation in x, y and z direction, 
respectively. It has an energy level spectrum with (25 + 1) = 21 values which, in first 
approximation, can be labelled by the quantum numbers M — —10, —9, ...10. The energy 
spectrum, can be obtained by using standard diagonalisation techniques of the [21x21] matrix 
describing the spin Hamiltonian S = 10. At H = 0, the levels M = ±10 have the lowest energy. 
When a field H z is applied, the energy levels with M « increase, while those with M » 
decrease. Therefore, different energy values can cross at certain fields. This crossing can be 
avoided by transverse terms containing S x or S v spin operators which split the levels. The spin 
S is in resonance between two states M and M' when the local longitudinal field is close to 
such an avoided energy level crossing (\H Z \ < 10~ 8 T for the avoided level crossing around H z 
= 0). The energy gap, the so-called tunnel spitting Am.m'j can be tuned by an applied field in 
the xy— plane via the S X H X and S y H y Zeeman terms. It turns out that a field in H x direction 
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(hard anisotropy direction) can periodically change the tunnel spitting A as displayed in Fig. [I] 
where if 2 in Eq. [j] was taken from Jl6| . In a semi-classical description, these oscillations are 
due to constructive or destructive interference of quantum spin phases of two tunnel paths . 

A direct way of measuring the tunnel splittings Am, At' is by using the Landau- Zener model 
| fl7[ |Tsfl which gives the tunnelling probability Pm,m> when sweeping the longitudinal field H z 
at a constant rate over an avoided energy level crossing Jl9| : 



P M M , = 1 — e 2 "»« l M - M ' l f dH » /<** (2) 

Here, M and M' are the quantum numbers of the avoided energy level crossing, dH z /dt is the 
constant field sweeping rate, g sa 2, the Bohr magneton, and h is Planck's constant. 

In order to apply the Landau-Zener formula (Eq. |^), we first cooled the sample from 5 K 
down to 0.04 K in a field of H z — -1.4 T yielding a negative saturated magnetisation state. 
Then, we swept the applied field at a constant rate over the zero field resonance transition 
and measured the fraction of molecules which reversed their spin. This procedure yields 
the tunnelling rate P_io,io and thus the tunnel splitting A_io,io (Eq. |^). The predicted 
Landau-Zener sweeping field dependence of P-io,io can be checked by plotting A_io,io as a 
function of the field sweeping rate which should show a constant which was indeed the case 
for sweeping rates between 1 and 0.001 T/s (fi g. B) . The deviations at lower sweeping rates 
are mainly due to the hole-digging mechanism |20|, [2l| which slows down the relaxation. The 
comparison with the isotopically substituted Fes samples shows a clear dependence of A_io,io 
on the hyperfine coupling (Fig. ^|). Such an effect has been predicted for a constant applied 
field by Tupitsyn et al. [|22| . 

All measurement so far were done in the pure quantum regime (T < 0.36 K) where transition 
via excited spin levels can be neglected. We discuss now the temperature region of small 
thermal activation (T < 0.7 K) where we should consider transition via excited spin levels 
[ p3| . We make the Ansatz that only ground state tunnelling (M = ±10) and transitions via 
the first excited spin levels (M = ±9) are relevant for temperatures slightly above 0.36 K. We 
will see that this Ansatz describes well our experimental data but, nevertheless, it would be 
important to work out a complete theory 

In order to measure the temperature dependence of the transition rate, we used the Landau- 
Zener method (l7) as described above with a phenomenological modification of the transition 
rate P (for a negative saturated magnetisation): 

P = n_i P_ 10,10 + Pth (3) 

where P_io,io is given by Eq. || , n_io is the Boltzmann population of the M = —10 spin 
level, and Pth is the overall transition rate via excited spin levels. n_io ~ 1 for the considered 
temperature T < 0.7 K and a negative saturated magnetisation of the sample. 

Fig. H displays the measured transition rate P for s *Fes as a function of a transverse field 
H x and for several temperatures. The oscillation of P are seen for all temperatures but the 
periods of oscillations decreases for increasing temperature (Fig. |]). This behaviour can be 
explained by the giant spin model (Eq. [l]) with higher order transverse terms (#2)- Indeed, 
the tunnel splittings of excited spin levels oscillate as a function of H x with decreasing periods 
(Fig-0). 

Fig. ^ displays the transition rate via excited spin levels Pth = P— "--lo-P-io.io- Surprisingly, 
the periods of Pth are temperature independent in the region T < 0.7 K. This suggests that 
only transitions via excited levels M — ±9 are important in this temperature regime. This 
statement is confirmed by the following estimation p6), see also Ref. p5|]. 
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Using Eq. [|, typical field sweeping rates of 0.1 T/s, and tunnel splittings from Fig. [l], one 
easily finds that the Landau Zener transition probability of excited levels are P-m,m ~ 1 
for M < 10 and H « 0. This means that the relaxation rates via excited levels are mainly 
governed by the lifetime of the excited levels and the time T res ,M during which these levels are 
in resonance. The later can be estimated by 

A_M,M 

Tres,M = 7T ,„ , , , ■ (4) 

gii B Mfj, dH z /dt 

The probability for a spin to pass into the excited level M can be estimated by T M 1 e~ El0 - M ^ kBl 
where E\q,m is the energy gap between the levels 10 and M, and tm is the lifetime of the 



excited level M. We yield (26 



P th W ^ T res,M g-fiip . M /k B T _ ^-M.M c -E 10 . M /k B T /qn 

Note that this estimation neglects higher excited levels with |M| < 8 |^7j. Fig. |^ displays 
the measured for the three isotopic Fes samples. For 0.4 K < T < 1 K we fitted Eq. [|to the 
data leaving only the level lifetimes rg and rs as adjustable parameters. All other parameters 
are calculated using Eq. |l] (28). We obtain t 9 = 1.0,0.5, and 0.3 x 10 _6 s, and r 8 = 0.7,0.5, 
and 0.4 x 10~ 7 s for ^Fes, s 'Fes, and 57 Fes, respectively. This result justifies our Ansatz of 
considering only the first excited level for 0.4 K < T< 0.7 K. Indeed, the second term of the 
summation in Eq. is negligible in this temperature interval. It is interesting to note that 
this finding is in contrast to hysteresis loop measurements on Mni2 [|30"f which suggested an 
abrupt transition between thermal assisted and pure quantum tunnelling fl3l|| . Furthermore, 
our result shows clearly the influence of nuclear spins which seem to decrease the level lifetimes 
tm, i.e. to increase dissipative effects. 

The nuclear magnetic moment and not the mass of the nuclei seems to have the major 
effect on the dynamics of the magnetization. In fact the mass is increased in both isotopically 
modified samples whereas the effect on the the relaxation rate is opposite. On the other hand 
ac susceptibility measurements at T > 1.5 K showed no clear difference between the three 



samples 32 suggesting that above this temperature, where the relaxation is predominately 
due to spin-phonon coupling |^3|, ^4|, the role of the nuclear spins is less important. Although 
the increased mass of the isotopes changes the spin-phonon coupling, this effect seems to be 
small. 

We can also exclude that the change of mass for the three isotopic samples has induced a 
significant change in the magnetic anisotropy of the clusters. In fact the measurements below 
T < 0.35 K, where spin-phonon coupling is negligible, have shown that (i) relative positions of 
the resonances as a function of the longitudinal field H z are unchanged |l3| , and (ii) all three 
samples have the same period of oscillation of A as a function of the transverse field H x || , a 
period which is very sensitive to any change of the anisotropy constants. 

In conclusion, we presented detailed measurements based the Landau Zener method which 
demonstrated again that molecular magnets offer an unique opportunity to explore the quan- 
tum dynamics of a large but finite spin. We believe that a more sophisticated theory is needed 
which describes the dephasing effects of the environment. 
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Fig. 1. - Calculated tunnel splitting A m ,m> (Eq. |l|) as a function of the transverse field H x for quantum 
transition between M = ±10, ±9 and ±8. 
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Fig. 2. - Field sweeping rate dependence of the tunnel splitting A_io,io measured by a Landau Zener 
method for three Fes samples, for H x = 0. The Landau Zener method works in the region of high 
sweeping rates where A_io,io is sweeping rate independent. Note that the differences of A_io,io 
between the three samples are rather small in comparison to the oscillations in Fig. j3|. 
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Fig. 3. - Transverse field dependence of the transition rate P at several temperatures, and the ground 
state transition rate P_io,io measured at T = 0.05 K and for s 'Fes. The field sweeping rate was 0.14 
T/s. The dotted lines indicate the minima of P_io,io- 
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Fig. 4. - Temperature dependence of the field between minima of the transition rate P in Fig. til a, 
b, and c are defined in the inset. The dotted line labelled with a', b', and c' where take from Pth of 
Fig. g; see also fl29| . 
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Fig. 5. - Transverse field dependence of P t u which is the difference between the measures tunnel 
probability P and the ground state tunnel probability n_io-P-io,io measured at T = 0.05K (see Fig. ^| 
The field sweeping rate was 0.14 T/s. The long dotted lines indicate the minima of Pth whereas the 
short dotted lines indicate the minima of P-io,io- 




a 57 F e 8 
o st Fe 8 
x D F e 8 

I 0.4K 



0.65K 



1.5 2 2.5 

1/T (1/K) 



Fig. 6. - Temperature dependences of Pth for H x = for three Fes samples. The field sweeping rate 
was 0.14 T/s. The dotted lines are fits of the data using Eq. | H|. 



